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ABSTRACT 

We present an analysis of the spatial distribution of star formation in a sample of 60 visu- 
ally identified galaxy merger candidates at z > 1. Our sample, drawn from the 3D-HST 
survey, is flux-limited and was selected to have high star formation rates based on fits of 
their broad-band, low spatial resolution spectral energy distributions. It includes plausible 
pre-merger (close pairs) and post-merger (single objects with tidal features) systems,with to- 
tal stellar masses and star formation rates derived from multi- wavelength photometry. Here 
we use near-infrared slitless spectra from 3D-HST which produce Ha or [Oiii] emission line 
maps as proxies for star-formation maps. This provides a first comprehensive high-resolution, 
empirical picture of where star formation occurred in galaxy mergers at the epoch of peak 
cosmic star formation rate. We find that detectable star formation can occur in one or both 
galaxy centres, or in tidal tails. The most common case (58%) is that star formation is largely 
concentrated in a single, compact region, coincident with the centre of (one of) the merger 
components. No correlations between star formation morphology and redshift, total stellar 
mass, or star formation rate are found. A restricted set of hydrodynamical merger simula- 
tions between similarly massive and gas-rich objects implies that star formation should be 
detectable in both merger components, when the gas fractions of the individual components 
are the same. This suggests that z ^ 1.5 mergers typically occur between galaxies whose gas 
fractions, masses, and/or star formation rates are distinctly different from one another. 

Key words: galaxies: formation - galaxies: interactions - galaxies: starburst - galaxies: struc- 
ture 



1 INTRODUCTION 

The spatial extent and distribution of star formation in normal, lo- 
cal galaxies is well established (e.g., James et al. 2004; Bigiel et al. 
2008, 201 1 ; Schruba et al. 201 1; Calzetti et al. 2012, and references 
therein). Both self-regulated star formation and merging must be 
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key ingredients in galaxy-formation and evolution and have been 
studied observationally in detail in the z < 0.5 Universe (e.g., 
Lambas et al. 2003; Hammer et al. 2005; Barton et al. 2007, 2000; 
Jogee et al. 2009; Robaina et al. 2009, 2010) as well as in theoreti- 
cal simulations (e.g., Mihos & Hernquist 1996; Barnes & Hernquist 
1996; Springel 2000; Cox 2004; Cox et al. 2006, 2008; Lotz et al. 
2008a,b; di Matteo et al. 2007). From these studies, it has become 
evident that rapid star formation can be triggered by tidal interac- 
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tion in mergers, but the simulations also suggest that mergers trig- 
ger both nuclear starbursts and black hole accretion. Even though 
galaxy mergers are observed to enhance star formation in galaxies 
and trigger some of the most violent starbursts known (e.g., Rieke 
et al. 1985; Joseph & Wright 1985; Melnick & Mirabel 1990; Klaas 
& Elsaesser 1991), it appears that the net effect of major mergers in 
the global star formation history of the galaxy population has been 
relatively modest since at least z = 1 (e.g., Robaina et al. 2009; 
Jogee et al. 2009). 

Mergers are predicted to play a crucial role in the build-up 
and formation of massive galaxies (e.g., Mihos & Hernquist 1996; 
Springel 2000; Cox et al. 2008; Hopkins et al. 2010), and there- 
fore, a crucial step towards fully understanding how galaxies have 
evolved is to study the star formation properties in merging systems 
at high redshifts. Some of the questions that need to be addressed 
observationally to get a more detailed picture of the star formation 
history of present-day galaxies are: Where did stars form in merg- 
ing galaxies at higher redshift? and Which phase(s) of the merging 
process seems to trigger most star formation? 

At higher redshift, the interplay between merging and star for- 
mation has been investigated in much less detail than in the lo- 
cal Universe. This is mostly due to observational challenges, as re- 
solved observations at < 1" resolution are needed. Further, tracing 
star formation, e.g., through Ha at z ^ 1.5 requires observations in 
the near-infrared (NIR), and these have been far less feasible than 
observations in optical wavebands. Nevertheless, the 0.7 < 2; < 2 
epoch is immensely important for understanding galaxy-formation, 
as this is the cosmic time when the majority of the stars we see in 
galaxies today were formed (e.g., Hopkins & Beacom 2006; Karim 
et al. 2011). Several studies have addressed the impact of (major) 
galaxy mergers as well as the general galaxy morphology on the 
amount of star formation at 2: ~ 1.5 (e.g., Swinbank et al. 2004; 
Law et al. 2007; Forster Schreiber et al. 2009; Wright et al. 2009; 
Conselice et al. 201 1; Block et al. 2012; V -V al. 2012). However, 
only a few of these studies investigate three-dimensional (3D) spec- 
troscopy, where both spatial and spectral information is available, 
which is crucial for investigating the spatial extent of star forma- 
tion. 

Until recently, large samples of galaxies, and in particular 
galaxy mergers, with rest-frame optical 3D spectroscopic informa- 
tion at high redshift did not exist. The largest samples of galax- 
ies with such data at 1.5 < z < 2.5 are from the SINS sur- 
vey (Forster Schreiber et al. 2009), recently expanded with the 
zC-SINF sample (Mancini et al. 2011), totalling 110 star-forming 
galaxies at 1.5 < z < 2.5 observed with SINFONI, and from 
the MASSIV survey (Contini et al. 2012) which contains 84 star- 
forming galaxies at 0.9 < z < 1.8 also mapped with SINFONI. 
Using Ha as kinematic and star formation tracer enabled analy- 
sis of the spatially-resolved ionised gas kinematics, its distribu- 
tion, and the physical properties of these systems (Genzel et al. 
2006, 2008, 2011; Shapiro et al. 2008; Cresci et al. 2009; Epinat 
et al. 2009, 2012; Queyrel et al. 2009, 2012; Newman et al. 2012). 
High-resolution NIR imaging with The Hubble Space Telescope 's 
(HST) Near Infrared Camera and Multi-Object Spectrometer 2 
(NICMOS2) for a small subset of the SINS objects provided addi- 
tional rest-frame optical morphologies, in agreement with the disc 
or merger nature from the Ha kinematics (Forster Schreiber et al. 
201 la,b). The selection of these samples was primarily based on in- 
tegrated photometry or spectroscopic properties, not morphologies, 
and only a modest fraction of objects (^-^1/3) were kinematically in- 
ferred to be (major) mergers. Roughly comparable fractions were 
found in other sizeable NIR IFU samples at 2: ~ 1 — 3, including. 



e.g., those by Law et al. (2007, 2009); Gnerucci et al. (2011), and 
Wisnioski et al. (2011b). 

With the recent 3D-HST slitless grism survey (see Brammer 
et al. 2012, and Section 2), much larger samples of objects with 
NIR 3D emission line spectroscopy have become available, mak- 
ing it possible to address the spatial extent of star formation for 
extensive samples of galaxy mergers at the peak of cosmic star for- 
mation rate density. The IFU samples mentioned above have signif- 
icantly higher spectral resolution, enabling detailed kinematic stud- 
ies, but AO-assisted IFU observations, which provide angular res- 
olution comparable to HST in the NIR, remain observationally ex- 
pensive and suffer from complications due to strong night sky lines. 
HST grism observations, as the ones taken in the 3D-HST survey, 
provide more limited kinematic information but allow for unbiased 
target selection and are much more efficient at detecting and map- 
ping the continuum and line-emission at high angular resolution for 
all targets within the field-of-view. 3D-HST provides resolved line- 
emission, enabling studies of the spatial extent of star formation 
for large samples of galaxies at 2; > 1 in five well- studied cosmo- 
logical fields. The initial papers use approximately half of the full 
dataset, as described in van Dokkum et al. (2011) and Brammer 
et al. (2012). 

Using the same 3D-HST data, we explore the spatial extent 
and distribution of star formation in a 'population snapshot' of pre- 
sumably merging systems at z ^ 1.5. For this sample we make 
~0.2^^ resolution maps of emission lines (Ha and [Olll]), which 
trace the spatial extent of the (unobscured) star formation in these 
mergers and allow us to study their star formation properties in a 
statistical and unbiased way. This is done under the assumption that 
the Ha (for z - 0.7 - 1.5) and [Oni] (for z - 1.2 - 2.3) emis- 
sion of the systems trace the star formation. This has been shown to 
be a fair assumption for both Ha (Kennicutt 1983; Gallagher et al. 
1984; Kennicutt et al. 1994; Kennicutt 1998a,b) and [Olll] (Kenni- 
cutt 1992; Teplitz et al. 2000; Hippelein et al. 2003), even though 
using [Olll] as a quantitative indicator of star formation rate (which 
is not what we aim to do here) includes several complicating factors 
(TepHtz et al. 2000). 

To help interpret our observations in a theoretical context, 
we create a sample of pseudo-observations from state-of-the-art 
smoothed-particle hydrodynamic (SPH) simulations of individual 
mergers, which we compare to the 3D-HST data. The simulations 
predict a merger sequence and star formation picture with centrally- 
concentrated triggered starbursts at final coalescence, enhanced star 
formation in tidal features, and black hole growth and accretion. 
The goal is to understand the observational results from 3D-HST 
by making a direct comparison with the predictions from the sim- 
ulations. These comparisons will help explore which parameters, 
e.g., viewing angle, merger phase, gas fraction, mass ratio, etc. play 
a crucial role in determining, for example, star formation rates from 
observations. 

In Section 2 we describe the 3D-HST survey from which the 
merger sample was selected. We then describe the selection of our 
sample in Section 3 and the procedure used to map the spatial ex- 
tent of the star formation in Section 4. In Section 5 we split the 
sample into four different morphological types and find that most 
mergers exhibit star formation in only one component. In Section 6 
we compare numerical merger simulations to the observed 3D-HST 
spectra and find that in simulations star formation most commonly 
occurs in both components, before we summarise and conclude in 
Section 7. 



© 2013 RAS, MNRAS 000, 1-16 



Spatial Extent ofSFinz ^ 1.5 Mergers 3 



2 THE 3D-HST SURVEY DATA 

To construct the 0.2" resolution emission line maps, the prox- 
ies for star formation maps, we take advantage of the NIR 3D 
spectroscopy survey possibilities that the Wide Field Camera 3 
(WFC3) on HST brings. The 3D-HST survey is a 248 orbit 
NIR spectroscopic Hubble treasury program (Cycles 18 and 19, 
PI van Dokkum). It provides NIR imaging with the F140W filter 
and grism spectroscopy with the G 141 grism over well-studied ex- 
tragalactic survey fields (AEGIS, COSMOS, GOODS-S, GOODS- 
N, and UKIDSS/UDS). The grism spectroscopy is slitless so both 
spatial and spectroscopic information is available for every single 
object in the survey fields. The 3D-HST survey provides rest-frame 
optical spectra for a sample of ~7000 galaxies at 1 < z < 3.5 (van 
Dokkum et al. 201 1; Brammer et al. 2012). 

As of August 2011 the survey had observed 68 pointings over 
the GOODS-S, G00DS-N\ COSMOS and AEGIS fields (van 
FioVViim q1 ^ni i)^ xhe present work is based on 30 of these 
68 pointings, where the extensive ancillary data available enables 
robust spectral energy distribution (SED) modelling necessary for 
our sample selection as described in Section 3. 



2.1 The 3D-HST Grism Spectroscopy 

The WFC3 G141 grism used in 3D-HST disperses the light over 
the wavelength range from 1.05 /im to 1.7 /im with a low spectral 
resolution ofR^ 130. 

Since the grism spectroscopy is slitless, the WFC3 G141 
grism basically produces an emission line image that is superim- 
posed onto a sequence of dispersed monochromatic continuum im- 
ages, and some of the key features of slitless spectroscopy there- 
fore need to be taken into account. First of all, the width of emis- 
sion/absorption lines in the dispersion direction in slitless spec- 
troscopy is not only caused by velocity broadening (which is neg- 
ligible for the low resolution 3D-HST spectra) and the intrinsic 
broadening of the wavelength dispersion: as slitless spectroscopy 
produces shifted monochromatic images, the spatial extent of the 
dispersed emission line image reflects the spatial distribution of the 
line emission both along and perpendicular to the dispersion di- 
rection. We will take advantage of this 'morphology-broadening' 
(which can be seen in the third panel of Figure 4), to map the spa- 
tial extent of star formation as described in Section 4. Figure 4 will 
be explained in more detail here. 

As with multi-slit spectroscopy, the differing wavelength cov- 
erage of the spectra is an issue. Since the detector onto which the 
field-of-view is dispersed has a finite size, approximately 10% of 
the spectra are cut off on the edge of the detector. 

Lastly, 'contamination' is an important property of shtless 
spectroscopy. Since the focal plane is not blocked out with a sht 
or a mask as is usually done in standard spectroscopy, all the light 
from a given object, and all other objects in the observed field, are 
dispersed onto the detector. Hence, spectra will often overlap and 
therefore 'contaminate' each other as explained in Brammer et al. 
(2012). 

For more information on the data reduction methods, the data 
products of the 3D-HST survey, and the survey itself, we refer to 
Brammer et al. (2012). 



1 The GOODS-N data were taken as part of the HST program GO- 1 1600 
(PI B. Weiner). 



3 SELECTING MERGER CANDIDATES 

We select our sample of merger candidates from the first 68 point- 
ings obtained as part of the 3D-HST survey based on three different 
inputs: (i) The 3D-HST survey catalogue, (ii) SED modelling, and 
most importantly (iii) visual inspection of NIR (F140W) morpholo- 
gies. The first two selection steps are performed to define an initial 
sample of systems with sufficient spectral coverage and to minimise 
the number of objects to visually inspect. We will describe each of 
these three steps below. 

To obtain robust star formation rate estimates via the SED fit- 
ting described in Section 3.2, we require extensive ancillary photo- 
metric catalogues. We therefore focus on the 30 pointings of data 
available in GOODS-S (6) and COSMOS (24), where the photo- 
metric data in the FIREWORKS (Wuyts et al. 2008) and the NEW- 
FIRM medium band survey (NMBS; Whitaker et al. 2011) cata- 
logues are available, respectively. Hence, this work is performed 
on approximately 1/5 of the final 3D-HST data product. 

3.1 Grism Catalogue Cuts 

The first step in defining our merger sample is to select a well- 
defined sample of objects based on the data products of the 3D-HST 
survey. We ensure that at least 75% of each spectrum in the sample 
falls on the detector. Since we are looking for merging objects we 
do not put any constraints on the contamination of the individual 
spectra, as spectra of close pairs will always have a high level of 
contamination. We rely on the visual inspection (Section ) to 
remove cases with heavy contamination from objects that are not 
part of the potentially merging system. 

Each individual object in the 3D-HST catalogue has been 
matched to the available ancillary photometric catalogues. Since 
the 3D-HST catalogue is selected from the deep (HmAow ~ 26.1; 
Brammer et al. 2012) high-resolution NIR HST F140W images, 
and the photometric catalogues are ground-based and shallower, 
not all 3D-HST objects can be matched to an object in the pho- 
tometric catalogues. We only selected objects with a counterpart 
(within 0.3'') in ground-based photometric catalogues. Faint ob- 
jects have the risk of being assigned to a bright(er) counterpart's 
photometry, have low signal-to-noise (S/N), and less reliable red- 
shifts, and we therefore restrict ourselves to objects with mpHow ^ 
23.5. 

Last but not least, the 3D-HST catalogues provide a redshift 
estimate for the objects based on the extracted grism spectra. The 
catalogue grism redshifts, Zgrism, are obtained by collapsing the 2D 
grism spectrum into a ID spectrum, combining it with available 
photometry, and then estimating the redshift with an updated ver- 
sion of the EAZY code (Brammer et al. 2008). The redshift range 
where Ha and/or [Olll] emission fall in the G141 grism wavelength 
range is 0.7 < z < 2.3 (see Figure 1 in Brammer et al. 2012). We 
are interested in tracing the star formation in the merging systems 
via either Ho; or [Olll] emission, and therefore use 2;grism to select 
objects in this particular redshift range. 

Applying these five initial cuts (listed in the top half of Ta- 
ble 1) reduces the full sample of 21460 detections in the 30 
GOODS-S and COSMOS pointings to 1542 objects. 

3.2 Fitting SEDs to Photometry 

We select star-forming systems that are expected to have significant 
emission line features based on their star formation rates (SFR), 
specific SFR (sSFR), and stellar masses (M*). We obtain the SFR, 
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Table 1. Grism (Top) and SED (Bottom) Selection Criteria 
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sSFR, and M* of each individual object from modelling the SED 
based on the ancillary photometric catalogues using a Chabrier 
(2003) IMF with the code MAGPHYS presented in da Cunha et al. 
(2008). We use the 37 NMBS bands (Whitaker et al. 201 1) for the 
COSMOS objects and the 17 FIREWORKS bands (Wuyts et al. 
2008) for the GOODS-S objects. Both catalogues span from the 
far-UV to MIPS 24 /xm. The 3D-HST catalogue redshift ^gnsm is 
used as a prior when fitting the ancillary photometric data for each 
object. 

Although the photometric measurements are a blend of two 
or more components in many cases, i.e., only one photometric ID 
corresponds to each merger whereas several components are clearly 
distinguishable in the high-resolution HST imaging, selecting the 
high-SFR objects based on SED fits to the photometry is still very 
effective in selecting objects with strong emission line features. 

As shown in the bottom part of Table 1 we select objects with 
SFR > lOMo/yr, sSFR > IQ-^ Vr"^ and lO^M© < M* < 
lO^^M©. A SFR of lOM© /yr roughly corresponds to an emission 
line flux Fucx 10~^^ erg/s/cm^ at z — 1.5 which corresponds 
to a (collapsed) emission line S/N of ~8 at 1 x 10~^^ erg/s/cm^ 
(Brammer et al. 2012). Hence, concentrated star formation, i.e., 
emitted from a modest amount of pixels of this order, should be 
well-detected in the 3D-HST spectra. On the other hand, if the to- 
tal emission line flux Fna is spread over a larger area, the S/N per 
pixel might become too low for clear detection (see Section 5). 

We find 352 of the initial 1542 objects in the 0.7 < z < 2.3 
range that satisfy these SED criteria. 

3.3 Visual Inspection of NIR Morphology 

Previous studies have argued both for visual (Robaina et al. 2009) 
and algorithmic merger identification (e.g., Lotz et al. 2008a; Con- 
selice et al. 2008). For this pilot study, which is focussed on the 
emission line morphology not on the merger rates, we have, as dis- 
cussed below, decided to use a visual classification. The visual in- 
spection of the remaining 352 objects is the crucial final step in the 
merger sample selection process. 

The visual inspection is based on the NIR F140W morphol- 
ogy of the objects from the 3D-HST direct imaging. The NIR im- 
ages show the rest-frame optical emission at the redshifts of our 
galaxies. We assume that the observed NIR (i.e., rest-frame opti- 
cal) morphology traces the distribution of the (intermediate-age) 
stellar component of the galaxies, and that it is therefore differ- 
ent from the stars being formed (current star formation) as traced 
by the emission lines. A caveat to such an assumption is that if a 
galaxy does not have a dominating intermediate-age stellar popula- 
tion and has a high SFR, the morphology in the rest-frame optical 
will to some extent reflect the distribution of young stars as well. 



The criteria used to select the merger candidates from the 352 ob- 
jects are that (i) they should show a morphology that differs from 
the bulk of the 'normal' isolated galaxies, i.e., a disturbed irregu- 
lar/asymmetric morphology, and (ii) they have to show several dis- 
tinct components in the continuum image, either multiple objects 
within the ~ 3 X 3 arcsec F140W thumbnails or pronounced tidal 
features extending from the main continuum emission component. 
It should be noted that because our merger selection is based only 
on this visual classification, and the low resolution of the G141 
grism {R ~ 130) does not offer any kinematic information of the 
individual companions, it is impossible to address whether or not 
the systems are gravitationally bound. Our merger sample therefore 
consists of potentially merging systems. 

A fraction of the more widely separated merger pairs could 
therefore be potential low- or high-redshift interlopers which would 
artificially enhance the number of mergers found. The distances be- 
tween the majority of the individual merger components in the can- 
didates selected here are of the order 1" , corresponding to roughly 
8.5 kpc dXz ^ 1.5. Law et al. (2012) estimates that lt\% of galax- 
ies have projected false pairs within 16 kpc. This serves as an up- 
per limit on the expected false pair fraction for our sample. As we 
also include late-stage mergers and not only widely separated pairs, 
based on this a more realistic estimate of the amount of interlopers 
would probably be ~ 4%. 

Furthermore, several studies have shown that the morphology 
of isolated star-forming galaxies at higher redshifts is often clumpy 
and irregular (Bournaud et al. 2007; Genzel et al. 2008, 2011; Law 
et al. 2009, 2012; Elmegreen et al. 2009; Kriek et al. 2009; Forster 
Schreiber et al. 201 la,b; Wuyts et al. 2012; Wisnioski et al. 201 la). 
This potentially biases our visual classification as inclusion of such 
systems will artificially enhance the number of selected mergers 
and hence the estimated fraction of mergers in our sample. As de- 
scribed below both parametric and visual classification schemes 
will be subject to this bias. Hence, what is described as a merger 
in the present selection might in fact be a galaxy with a clumpy and 
irregular light distribution appearing like a merger remnant. This 
caveat should be kept in mind when evaluating the merger candi- 
dates presented here and elsewhere. In the remainder of the paper 
we will therefore use the shorthand merger for likely merger can- 
didates. 

The visual inspection is not only used to select morphologi- 
cally disturbed systems. Inspecting the full grism spectra, as well 
as the one-dimensional collapsed grism spectra, of the individual 
objects, we are able to discard objects without any emission line 
features. Without emission line features creating a star formation 
map as described in Section 4 is impossible. Assuming that the 
estimated SFR, mass, and redshift are correct, the fact that these 
objects lack emission line features make them very interesting in 
themselves, as they might be highly dust-obscured systems block- 
ing the star formation emission. However, these objects are ignored 
for the present study. If an emission line feature on the other hand 
is observed, i.e., the (collapsed) emission line flux is roughly larger 
than ~ 10~^^ erg/s/cm^ (see Section 3.2), the object is included in 
our sample and we attempt to create an emission line map. How- 
ever, as noted in Section 3.2 and as we will see in Section 5 this 
does not necessarily mean that the S/N per pixel in the full grism 
spectrum is good enough to create an emission line map. 

Visually inspecting the grism spectra also ensures that any 
strong contamination is due to the different components of the 
merging system and not due to interlopers. 

From the parent sample of 352 catalogue-selected objects, the 
visual inspection discards 292 objects as they appear to be iso- 
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lated 'normal' galaxies (252/292), or having high contamination 
not stemming from the merging components (13/292), or having no 
obvious emission line features in their spectra (24/292 correspond- 
ing to 7^2% of the 352 catalogue- selected objects). Hence, our fi- 
nal sample of (potential) mergers from 6 GOODS-S and 24 COS- 
MOS 3D-HST pointings consists of 60 systems, corresponding to 
^ 17^4% of the 352 catalogue- selected objects, or approximately 
10-13% if we correct for the expected fraction of false pair inter- 
lopes. Here the confidence intervals are the 95% quantiles of the 
beta distribution following Cameron (201 1). This merger candidate 
fraction is on the high side compared to what is generally found in 
the literature (see, e.g.. Figure 1 in Lotz et al. 2011 and Williams 
et al. 2011), again suggesting that a fraction of the visually iden- 
tified mergers may be single objects with a blotchy distribution of 
young stars which increases the apparent merger fraction. The ob- 
tained merger fraction is subjective as it relies on a visual assess- 
ment of disturbance. As described above, also potential galaxy in- 
terlopers which have passed the subjective visual inspection could 
be biasing our sample towards higher merger fractions. Williams 
et al. (2011) used a mass-selected sample (log(M/M0) > 10.5) 
of galaxy pairs when estimating a major merger fraction of ^5%. 
This could also account for part of the discrepancy, as we are not 
only looking at distinct pairs of galaxies and use a flux limit as 
opposed to a mass limit in the selection process. With an average 
merger fraction just below 10% at ^ ^ 1 the merger fractions sum- 
marised in Figure 1 of Lotz et al. (2011) which includes samples 
selected on both mass and luminosity cuts is somewhat closer to 
the 17^4% (^10-13% if corrected for interlopes) reported here. 

The selection cuts in Table 1 positions all 60 3D-HST mergers 
in the blue cloud of star-forming galaxies (Bell et al. 2004; Strateva 
et al. 2001) in colour magnitude diagrams, which implies (as ex- 
pected) that the selected mergers are mainly gas-rich 'wet' mergers. 
We note however, that this may be very different for mass-selected 
samples where no limits have been imposed on the overall SFR of 
the objects. 

In Figure 1 we show the mpuow magnitude distribution of the 
60 3D-HST merger candidates. In Figure 2 and 3 we plot them to- 
gether with the selection regions from Table 1 (grey shaded regions) 
as large solid symbols. The small grey points represent the 292 ob- 
jects discarded by the visual classification, i.e., the general galaxy 
populations satisfying the selection cuts in Table 1 . In Figure 3 the 
distributions of SFR, sSFR, %ism, and M* for the 60 merger can- 
didates are shown as histograms on the axes of the scatter plots. 
In both the histograms in Figure 1 and in Figure 3, the dotted lines 
correspond to the 16th and 84th percentiles of the distributions, and 
the dashed lines show the median values. In Section 6 we will use 
these values to determine the parameter space to sample when sim- 
ulating 3D-HST grism spectra. 

Visual inspection is a more subjective way of selecting merg- 
ers than for instance empirically established parametric merger 
classification schemes like for instance the Gini (C), M20 and 
CAS selections (e.g., Lotz et al. 2004, 2006, 2008a,b; Conselice 
2003; Conselice et al. 2008, 2009; Papovich et al. 2005; Scarlata 
et al. 2007). Both parametric and visual classification schemes have 
advantages and disadvantages. Determining a merger population 
based on parameters ensures that the selection is done in a consis- 
tent and uniform way for all objects. However, at redshifts where 
the NIR morphology is not fully understood and where galaxies 
look clumpy and irregular as described above, and where disturbed 
morphologies are prominent without necessarily being part of a re- 
cent merger, a parametric classification scheme might fall short of 
a visual classification. On the other hand visual classification is po- 
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Figure 1. The cumulative distribution of mpuow magnitudes for the 60 
3D-HST mergers. The dotted lines indicate the 16th and 84th percentiles of 
the distribution whereas the dashed line shows the distribution median. The 
16th and 84th percentile values are used when defining the parameter space 
for the simulated grism spectra in Section 6. 
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Figure 2. The distribution of inferred stellar mass as a function of redshift 
for the 60 3D-HST merger candidates (large points). The small grey points 
represent the 292 objects discarded by the visual inspection and the grey 
shaded region shows the selection region from Table 1 . 



tentially biased by the subjectivity of the classifier. Nevertheless, 
the human eye is known to be excellent at detecting and distinguish- 
ing features in noisy images and spectra and arguably minimises the 
bias of the clumpy irregular morphology of high- 2; systems. 

To ease comparison with the extensive literature using para- 
metric merger classifications we have estimated the G, M20, C, 
A, and S morphological parameters for the 60 visually selected 
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Figure 3. The 60 3D-HST mergers (large symbols) plotted in the main region of the catalogue selection space given in Table 1 (grey shaded regions). The 
small grey points represent the 292 objects discarded by the visual inspection. The hashed regions are empty due to our selection criteria. The histograms 
attached to the scatter plots show the distribution of Zgrism, M^, SFR, and sSFR for the 60 mergers. The dotted and dashed lines in the histograms indicate 
the 16th and 84th percentiles and the median of the distributions, respectively. These values are used to define the parameter space of the simulated grism 
spectra described in Section 6. The different symbols (and colours) indicate the morphology of the emission line maps (SF type) described in Section 5 as 
indicated above the panels. No obvious trends are found between the SF morphology and the SFR, sSFR, redshift, or stellar mass. 



merger candidates and the 292 visually discarded objects as shown 
in Appendix A. The selected candidates partially satisfy the empir- 
ical parametric merger selection but in general seems to be an av- 
erage subset of the parent distribution, i.e., not distinguishing itself 
clearly from it. Assuming that the 60 merger candidates are reli- 
able therefore speaks in favour of using visual classification when 
the NIR morphology is complicated, as a parametric GM20CAS 
classification would not be able to clearly distinguish the 60 merger 
candidates from the parent population. 

Despite the limitations of a visual selection of mergers and the 
only partial agreement with the empirical parametric classifications 
(Appendix A), we believe that the advantages of the visual classifi- 



cation scheme described above outweighs a 'blind' parametric se- 
lection for a study like the one presented here. Furthermore, we are 
probing an unexplored regime (NIR at high redshift) and the cur- 
rent parametric methods might not be appropriately calibrated and 
tested here and we have therefore chosen to use the visual selection 
in the reminder of this paper. 



4 EMISSION LINE MAPPING 

To quantify the extent of (unobscured) star formation in the 60 3D- 
HST mergers described in the previous section, we rely on the spa- 
tial information of the Ha and [Olll] emission lines that the slit- 
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Figure 4. Steps performed to obtain the emission line (star formation) maps 
as described in Section 4. The top panel shows a standard 3D-HST spectrum 
with a prominent emission line feature. The 2D continuum model, a set of 
polynomial fits to the observed spectrum when excluding the emission line 
feature (Section 4.1), is shown just below that. The third panel shows the 
spectrum from the top panel after subtraction of the 2D continuum model. 
The red square marks the emission line map cut-out at Accmax obtained 
by cross-correlating the NIR F140W thumbnail of the object (grey scale 
bottom panel) with the continuum-subtracted spectrum. The bottom panel 
illustrates how this emission line map cut-out is mapped back (red contours) 
onto the continuum image. 



less grism spectroscopy provides. From the grism spectra we cre- 
ate emission line maps by subtracting a model of the continuum 
light in the grism spectra such that only the probed emission line 
feature is remaining. This can then be mapped back onto the NIR 
continuum light-distribution of the object. In practice we: 

(i) Create a 2D continuum model for the grism spectrum. 

(ii) Subtract this continuum model from the 2D grism spectrum 
itself. 

(iii) Cross-correlate the F140W thumbnail with the continuum 
subtracted spectrum to map the (cut-out) continuum- subtracted 
thumbnail back onto the F140W thumbnail. 

Each of these steps is described in detail in the following subsec- 
tions and are illustrated in Figure 4. 



4.1 Continuum Modelling and Subtraction 

The 2D continuum models we subtract from the 3D-HST spec- 
tra are based on a third order polynomial fit to a one-dimensional 
spectrum. The ID spectra are obtained from a weighted sum of 
the individual lines in the 2D spectra. The polynomial fit to the 
ID spectrum is turned back into a 2D continuum model by con- 
catenating rows with the ID polynomial form weighted by a 'slit- 
profile' obtained from the columns blue- ward and red- ward of the 
probed emission line feature in the full 2D spectrum. Subtracting 
this model from the 3D-HST spectrum returns a two-dimensional 
emission line map as illustrated in Figure 4, where all that is left is 
the emission line feature. This approach is similar to the one used in 
Nelson et al. (2012). Often when dealing with slitless spectroscopy, 
and the 3D-HST grism spectra in particular, the goal is to remove 
contamination in a systematic manner. However, we relied on vi- 
sual inspection to remove badly contaminated objects, since merg- 
ers per definition are contaminated. Remaining contamination not 
affecting the continuum flux of the merger was masked out when 
modelling the continuum before subtraction. 



4.2 Constructing Emission Line Maps 

We used a simple cross-correlation between the NIR F140W 
thumbnail continuum image of each object and the corresponding 
full 2D emission line map in order to map the emission line map 
back onto the NIR image. In practice we calculate 



- ^widt h - ^widt h 

CC(A) = /i,j,F140W X /i,j,ELmap 



(1) 



for each of the first k = A^2d — A^width columns in the full 2D 
emission line map, where A^2d is the number of columns in the 2D 
emission line map and A^width is the width of the F140W thumbnail. 
The /z,j,Fi4ow and /i,j,ELmap is the flux in the pixel (i, j) for the 
F140W thumbnail and 2D emission line map cut-out (indicated by 
the red box in Figure 4), respectively. 

The maximum of the cross-correlation function, CC, indicates 
the wavelength, Accmax, where there is the largest overlap between 
the NIR light distribution of the object and the /cth cut-out of the 
full emission line map. The Accmax corresponds to the /cth column 
of the 2D emission line map plus A^width/2. 

In Figure 5 we show a collection of emission line maps (red 
contours) from our 3D-HST merger sample. The individual maps 
correspond to the region at Accmax that has been mapped back onto 
the NIR continuum image as illustrated in Figure 4. In each map 
we have marked the relative location of [Nil] AA6548,6583, Ha 
A6563, [Sii] AA6716,6730 andH/3 A4861, [Oiii] AA4959,5007 for 
the Ha and [Olll] maps, respectively. In some cases the [Olll] dou- 
blet is marginally resolved as seen in the upper left emission line 
map in Figure 5. The Ha- [Nil] composite is however not resolved 
in the 3D-HST grism resolution. We do not attempt to de-convolve 
these emission lines when creating the emission line maps. Assum- 
ing that F[oiii] A5007 = 3 X F[oiii] A4959 the redshift uncertainty im- 
posed by ignoring the [Olll] doublet is only 0.0024. This is less 
than the quoted average 3D-HST redshift precision of 0.0034(l+z) 
(Brammer et al. 2012) and is therefore not affecting the conclusions 
of this study. 

The different 'morphologies' of the emission line maps in Fig- 
ure 5 will be addressed in Section 5. In Appendix B we show the 
full sample of 3D-HST merger emission line maps. 
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Figure 5. Examples of Ho; and [Oni] emission line maps (red contours) plotted on the WFC3 F140W thumbnails (grey scale) for a sub-subsample of the 60 
3D-HST mergers. The dashed boxes indicate the mapped region. Each row represents one morphological type of star formation distribution as described in 
Section 6.2. From top to bottom each row show maps of SF type 1 ('one component' maps), SF type 2 ('both (all) component' maps), SF type 3 ('in-between' 
maps), and SF type 4 ('low S/N' maps). The Type 1 and 2 maps (first two rows) have been divided into into pre- (left) and post-mergers (right). The vertical 
white lines indicated the relative distance between [Nn] AA6548,6583, Ua A6563, [Sn] AA6716,6730 and H/5 A4861, [Oni] AA4959,5007 for the Ua and 
[Oni] maps, respectively. The full sample of emission line maps is shown in Appendix B. 



5 RESULTS: THE SPATIAL EXTENT OF STAR 
FORMATION IN HIGH-Z MERGERS 

As noted in the introduction, star formation and mergers are im- 
portant parts of understanding how high-z galaxies evolved into 
the galaxies we observe in the low-^ Universe. In the previous sec- 
tions we have described how we select the mergers from the 3D- 
HST data, and subtract the continuum light in the spectra to create 
the emission line (star formation) maps. In this section we charac- 
terise the morphological type of the spatial extent of the star forma- 



tion in the 60 3D-HST mergers. However, we first divide our sam- 
ple into two different kinds of mergers: the 'pre-mergers' and the 
'post-mergers'. By pre-mergers we mean objects that show mul- 
tiple clearly-distinct and pronounced continuum peaks in the NIR 
images, i.e., the optical continuum emission comes from multiple 
objects in the process of merging or about to merge. Examples of 
those are shown to the left in the two top panels of Figure 5. The 
post-mergers, on the other hand, are systems that have (presum- 
ably) undergone merging and now appear to be dominated by a nu- 
clear feature in the continuum with pronounced tidal features sur- 
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rounding it. Examples of these systems are shown to the right in the 
two top panels in Figure 5. Dividing the 60 mergers into these two 
sub-samples return 32 pre-mergers and 28 post-mergers. We note 
that the distinction between pre- and post-mergers is (operatively) 
our distinction to guide-the-eye when classifying and inspecting 
the emission line maps of the individual objects as described be- 
low. The identification with the merger phases is plausible, but will 
need more modelling. 

To characterise the location and the spatial extent of the star 
formation in the 3D-HST mergers we categorise the star formation 
maps into the following four morphological types (SF type): 

1) One Component: The star formation in the primary emis- 
sion line feature is significantly stronger than any secondary emis- 
sion line feature. The threshold used is Fp > 2.5 x F^, where Fp 
and Fs is the estimated aperture flux of the primary and secondary 
emission line feature, respectively. For the pre-merging systems the 
primary emission line feature corresponds to one of the multiple 
objects and for the post-mergers it refers to either the nuclear re- 
gion or a tidal feature. 

2) Both (All) Components: The star formation is pro- 
nounced/detected in all (or the majority if more than two) com- 
ponents of the system, i.e., Fp < 2.5 x Fs. 

3) In-between: The mapped star formation appears to be 
emerging from in-between the merging components. None of the 
post-mergers show this feature, so indeed this means in-between 
clearly distinguishable objects. 

4) Low S/N: The S/N per pixel of the emission line features in 
the 2D grism spectrum is too low to produce a convincing emis- 
sion line map. This may be the case for extended star formation 
as pointed out in Section 3.2. It can also be due to very dust- 
enshrouded star formation making the emission lines very weak. 

Each of the four rows of star formation maps in Figure 5 show 
examples of these four SF types. The results from characterising 
the two classes of mergers with these four SF types are shown in 
Figure 6. The error-bars are obtained by bootstrapping the results, 
i.e., by randomly drawing 60 SF types from the results 1000 times 
and then using the 2a width of the resulting SF type distributions 
as error-bars (hence no error-bar on the post-merger SF type 3 in 
Figure 6, as none were found). For both the pre- and post-mergers, 
the star formation is most prominent in just one of the components 
(SF type 1) for roughly 3/5 of the objects. In roughly 1/3 of the 
objects, star formation was detected in all components (SF type 
2). Hence, the distribution of the spatial extent of star formation 
among the pre- and post-mergers is consistent. In Table 2 we have 
listed the fractions for all 60 3D-HST mergers resulting from the 
classification of the emission line maps. 

The difference in the rates of objects with prominent star for- 
mation in just one component (SF type 1) and and mergers with 
star formation of type 2 might be a consequence of dust obscura- 
tion. We are only able to probe the unobscured star formation, so 
in cases where one component (or the tidal feature) is much more 
dust-obscured than the other we would end up with star formation 
maps of type 1. The discrepancy between the rates could also be 
due to different SFRs in the different components. Since the merg- 
ers are selected based on morphology and we do not have any kine- 
matic information, the fraction of SF type 1 objects might be biased 
by chance superpositions of objects on the sky at different redshifts, 
such that we only see line emission from one object in the NIR. As 
described in Section 3.3 we expect approximately 4-7% of such 
objects due to the modest distances of ~8.5 kpc involved. 

The results indicate that most mergers happen between objects 
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Figure 6. The morphological types of the star formation distribution (Sec- 
tion 5) observed in the 60 3D-HST merger candidates split into 'pre- 
mergers' of multiple individual objects (32 objects, black stars) and 'post- 
mergers' of systems with a nucleus and tidal features (28 objects, red 
squares). The number of objects with a given SF type is indicated above 
each point. The error-bars are obtained via bootstrapping as described in 
the text. The two samples have similar star formation distributions. About 
20% of the objects show star formation in all merger components (individ- 
ual objects or nucleus and tidal feature) whereas ~60% of the systems only 
show star formation in one component. 



Table 2. The Spatial Extent of Star Formation 



SF type 



3D-HST 



Simulations 



1) One Comp. 


58^^^% 


(35/60) 


28l^% 


(83/296) 


2) Both (All) Comp. 


32l^^% 


(19/60) 


59+^% 


(175/296) 


3) In-Between Comp. 


3l^% 


(2/60) 


0+g% 


(0/296) 


4) Low S/N per pixel 




(4/60) 


13l^% 


(38/296) 



Uncertainties are obtained via bootstrapping as described in the text. 
See Figure 8 for a plot of these values. 



of different gas fractions and/or different SFR, i.e., two merging 
components with significantly different properties. We will show 
below that this is backed up by an initial comparison with simulated 
mergers. 

In Figure 3 the different SF types are represented by different 
symbols to look for dependencies between the morphology of the 
star formation maps and SFR, sSFR, %ism, and M*. These quan- 
tities are obtained on the total photometry of the merging compo- 
nents as only a single photometric ID was assigned to each merger 
and hence includes the flux of the total system. The absence of 
correlations seen in Figure 3 suggests that all star formation mor- 
phologies occur at all redshifts irrespective of SFR and mass. 
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6 SIMULATING 3D-HST SPECTRA 

With the exceptional data of the 3D-HST merger sample pre- 
sented above, we can perform comparisons with the star forma- 
tion produced in simulated mergers at high redshift. Current high- 
resolution merger simulations are able to predict the spatial distri- 
bution of star formation in mergers and produce simulated images 
by including the emission by young, newly-formed stars and the 
transfer of starlight through gas and dust (Cox et al. 2006, 2008; 
Jonsson et al. 2010). Lotz et al. (2008a,b, 2010a,b) used these and 
similar simulations to describe the correlations between the rest- 
frame optical morphology of mergers, their mass ratios, total star 
formation, projected size, gas fractions, and merger time-scales, en- 
abling comparisons with observations. 

Simulations have shown examples of star formation triggered 
by direct galaxy interaction, star formation originating in the cen- 
tral cores of the individual merging components, and star formation 
appearing in tidal features. In the nearby Universe, many examples 
of such features have been found, e.g., in the Antennae Galaxies 
(Wang et al. 2004). Whether the predictions of simulations are also 
representative of the star formation in mergers at 2: ~ 1.5, i.e., 
when the Universe was only 3-4 Gyr old, have not been tested 
yet. With 3D-HST, the sample size of data with both rest-frame 
optical morphology and star formation morphology is becoming 
large enough that we can start looking at a population (snapshot) of 
mergers instead of individual case studies, and hence comparisons 
with predictions from simulations become feasible. In this section 
we make an initial illustrative attempt at comparing state-of-the-art 
merger simulations with the observational results from 3D-HST. In 
the following we describe how we create simulated WFC3 G141 
grism spectra from a small initial set of high resolution simulated 
mergers. 



6.1 SPH Merger Simulations 

Knowledge about both the spatial and spectral extent of the objects 
is crucial for simulating grism spectra. As input for our grism simu- 
lations we use three-dimensional data cubes of simulated mergers. 
These input data cubes are based on the simulations described in 
Younger et al. (2008) and are designed in the spirit of the Cox et al. 
(2006, 2008) simulations. The simulations are AT-body/SPH simu- 
lations from the GADGET code (Springel et al. 2001) from which 
simulated observations are generated with the radiative transfer 
code SUNRISE (Jonsson 2006; Jonsson et al. 2010). For our ini- 
tial comparison of simulations with actual data, we use the three 
models listed in Table 3. The three simulations are all of mass ra- 
tio 2:1 and the two merging objects in each simulation both have a 
gas fraction (/gas) of 40%, making them 'wet' mergers like the 3D- 
HST systems and mimicking the suggested gas fraction of 2: > 1 
galaxies (Daddi et al. 2010; Tacconi et al. 2010). The three mergers 
happen at three different relative orbital inclinations (G) of the two 
merging components. Hence, what we test here is the effect of ori- 
entation on the detected star formation distribution. We further have 
a series of time steps or snapshots of each galaxy (Nt) and a set of 
viewing angles (A^view) that we can simulate spectra for. Thus, the 
simulation parameter space sampled here is spanned by time, orien- 
tation of the two mergers, and viewing angle. As described below, 
we combine this with a set of parameters determined by the 3D- 
HST sample to define the total parameter space to be covered by 
the simulated spectra. 



Table 3. The Simulated 3D Input Data Cubes 



Name 


Nt 


N ■ 


Mass ratio* 


/gas 


e 


Sim. 1 


6 


3 


(Sb) 2:1 (Sc) 


0.4:0.4 


30° 


Sim. 2 


5 


3 


(Sb) 2:1 (Sc) 


0.4:0.4 


90° 


Sim. 3 


3 


3 


(Sb) 2:1 (Sc) 


0.4:0.4 


150° 



* Parenthesis give disc model (Younger et al. 2008). 



6.2 Simulating 3D-HST grism Spectra 

Simulating the grism spectra from the merger simulation output 
data cubes presented above is fairly straightforward. The mock 
grism spectra are created by dividing the data cube into wavelength 
slices, i.e., images corresponding to a 'filter' of width AA. Offset- 
ting or dispersing this sequence of images and co-adding the fluxes 
results in a grism spectrum. Using this approach, we turn the in- 
put data cubes into a sequence of simulated WFC3 G141 grism 
spectra. We use a pixel scale of 0.06^' and a spectral resolution of 
AA = 22A according to the 3D-HST grism spectra characteris- 
tics described in Brammer et al. (2012). For the HST PSF we use 
a Tiny Tim PSF^ and for the system throughput we used the G141 
sensitivity curve. 

We use the parameter space of the observed merger sample 
from 3D-HST to define the parameter space of the simulations. By 
re- scaling each data cube we simulate grism spectra corresponding 
to the 16th and 84th percentiles of the magnitude (mpwow), redshift 
(zgrism), and SFR distributions of the 3D-HST data as indicated by 
the dashed lines in the corresponding histograms in Figure 1 and 3. 

To adjust the SFR of the data cubes before turning them into 
grism spectra, we assume that the SFR to Ha-flux conversion fol- 
lows the empirical relation of Kennicutt et al. (1994): 

SFR ^ 1 Lhc 
[Mo/yr] 1.26 x lO^i [erg/s] ' ^ ^ 

This relation assumes a Salpeter (1955) initial mass function and 
is applied to the emerging SFR, i.e., after the effects of extinc- 
tion are taken into account by SUNRISE. The ttifhow is obtained 
by scaling the flux (including emission lines) integrated over the 
F140W passband independently. Combining the values from Ta- 
ble 3 with these data-determined parameters, we end up with a 
sample of 336 simulated spectra spanning the parameter space 
(Nt, A^view, B, mpHow, z, SFR). 

Added to the raw mock grism spectra are Poisson noise and 
noise terms corresponding to the read noise and dark current quoted 
on the WFC3 website, as well as the background sky levels pre- 
sented in Brammer et al. (2012). 

In Figure 7 we present a sequence of simulated G141 grism 
spectra. The spectra shown have various combinations of Nt, A^view, 
6, mFi4ow, and SFR. The thumbnails on the left are the noise- 
free simulated F140W images of the objects. 

6.3 Creating Star Formation Maps of Simulated Spectra 

To create the star formation maps for the simulated spectra, we 
treat them in exactly the same way as the actual 3D-HST spec- 
tra (see Section 4). First, we visually inspect the spectra to make 
sure that an emission line feature is available; for some of the faint 
magnitude, low-SFR combinations the emission lines do not show 

^ http://www. stsci.edu/hst/observatory/focus/TinyTim 
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Figure 7. Examples of simulated 3D-HST data based on simulated observations of N-body/SPH mergers generated with the radiative transfer code SUNRISE. 
The thumbnails (left) are emission line maps similar to the ones shown for the 3D-HST data in Figure 5. The first three spectra (right) result in star formation 
maps of SF type 1, while the last three result in SF type 2 star formation maps. Various combinations of SFR, mpHow. time-step in the merging of the 
components {Nt), viewing angles (A^view). and initial relative orbital inclination of the two merging components (B) are shown. 
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up in the noise-added spectra. Having eliminated spectra without 
emission Une features, as well as cases of viewing angle and time 
step were the two merging objects could not be distinguished, the 
sample of 336 simulated spectra is reduced to 296. For these 296 
spectra, we create emission line maps by the method described in 
Section 4. 

As for the actual 3D-HST data, we characterise the 296 emis- 
sion line maps of the simulated spectra according to the four SF 
types described in Section 5. The results from this classification are 
shown in Table 2 together with the results for the 3D-HST mergers. 
The simulated spectra in Figure 7 show examples of both spectra 
resulting in star formation maps of type 1 (first three spectra) and 
type 2 (last three spectra). 

6.4 Comparing the 3D-HST Data with Simulations 

We compare the 3D-HST SF types we inferred in Section 5 with the 
simulated SF types in order to investigate how, e.g., SFR and view- 
ing angle influences the frequencies of different types of mergers 
from the characterisation of the extent of star formation in mergers. 
In Figure 8 the SF types of the full sample of 60 3D-HST mergers 
(blue circles) are shown together with the SF types of the 296 sim- 
ulated spectra of the 3 mergers from Table 3 (green triangles). The 
plotted percentages are given in Table 2. The error-bars are again 
obtained via bootstrapping. The comparison needs to be done with 
care, as the sample of simulated mergers from Table 3 does not 
span the gas fraction and mass ratio dimensions of the parameter 
space, expected to influence the 'observability' of star formation 
(Lotz et al. 2010a,b, 2011), as they are fixed at 40% and 2:1, re- 
spectively (see Section 6.1). Nevertheless, our analysis shows how 
the 3D-HST survey enables a direct comparison of predictions from 
simulations with an extensive sample of high-z mergers. 

Despite the limitations of he spanned simulation space, it 
seems that the fraction of spectra where the S/N per pixel is too low 
to create an actual emission line map is fairly consistent. Likewise, 
the fraction of cases where the star formation seems to emerge from 
in-between the merging components (SF type 3) is comparable: in 
fact we do not find any objects in this category among the simu- 
lated spectra. Looking at the two SF type 3 cases we find for the 
data (third row of Figure 5) could indicate that the reason we are 
not seeing star formation from the components themselves is due 
to dust-obscuration rather than more pronounced star formation in- 
between the components. To find similar trends in the simulated 
emission line maps we would need to investigate a range of gas 
fractions rather than just the /gas = 40% cases presented here. This 
is however beyond the scope of this initial comparison. 

From Figure 8 it is also evident that the simulated spec- 
tra, on average, have twice as many star formation maps where 
both components have pronounced star formation (SF type 2), i.e., 
nearly 2/3 of the simulated maps, as compared to ^1/3 for the 3D- 
HST star formation maps. Correspondingly, the fraction of single- 
component emission line maps is significantly lower for the simu- 
lations than for the observations. As the simulations assume fairly 
similar properties between the two merging components (same gas 
fraction, similar mass in the form of major mergers as opposed to 
minor mergers, same dust content, and total SFR), this suggests 
that the majority of mergers do not happen between such galaxies. 
According to the 3D-HST data, only ^1/3 of real mergers seem to 
happen between galaxies of similar properties. For example, differ- 
ent gas fractions or a larger gap between the masses (minor merg- 
ers) could change this picture, in agreement with the results pre- 
sented in Section 5. 
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Figure 8. Comparing the spatial extent of star formation (SF type; Sec- 
tion 5) in the 60 3D-HST mergers (blue circles) with the 296 simulated 
spectra (green triangles). Each point has been assigned an error-bar obtained 
via bootstrapping. Direct comparisons should be done with care due to the 
limited size of the parameter space the simulations span. Nevertheless, it is 
evident that the simulations show many more cases where star formation is 
seen in both the merging components than the 3D-HST data do, indicating 
that the majority of mergers have different gas fractions and/or mass ratios 
prior to merging. 



Also, the distribution of orbits and orientations of the merger 
components could play a significant role here. The presented simu- 
lations consist of both prograde-prograde (Sim. 1), prograde-polar 
(Sim. 2) and prograde-retrograde (Sim. 3) mergers. We find a deficit 
of SF type 1 mergers among the Sim. 1 emission line maps, in good 
agreement with the notion that prograde-prograde mergers gener- 
ally result in more symmetric star formation (di Matteo et al. 2007). 
Hence, several factors might play a role in the difference found be- 
tween the 3D-HST and simulated emission line (star formation) 
maps. 



7 CONCLUSION 

We have presented a sample of galaxy merger candidates with 
three-dimensional (R.A., Dec, and A) spectroscopy at redshift 
z ^ 1.5. The sample consists of 60 morphologically selected merg- 
ers from the Hubble treasury slitless grism survey 3D-HST with 
total masses and star formation rates for the systems derived from 
multi- wavelength photometry. From the slitless grism spectroscopy 
we created emission line maps of the rest-frame optical emission 
lines Ha and/or [Olll] as a proxy for the spatial extent of the un- 
obscured star formation. Our results go towards a comprehensive 
empirical picture of where star formation happens in galaxy merg- 
ers at z ^ 1.5 when the cosmic star formation and merger rate were 
at their peak. 

We have also carried out an illustrative comparison of the 3D- 
HST mergers to recent SPH simulations of galaxy mergers which 
include star formation and dust extinction. This small sample of 
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simulations points the way toward more extensive simulation pro- 
grams aimed at spanning the full parameter space of the observa- 
tions. 

The main conclusions of the present study are: 

• The spatial distribution of star formation in z ^ 1.5 candidate 
mergers shows a broad range of morphologies. It is often concen- 
trated in a single, compact region, but can also be located in tidal 
tails or in-between the main stellar bodies of the progenitors. 

• In the majority (581^3%) of the early-stage, pre-coalescence 
mergers, the star formation is significantly more pronounced in just 
one of the merging components. This is likely due to different gas 
fractions of the progenitors. Alternatively, dust content may differ 
in quantity or configuration. 

• The star formation morphologies show no clear correlations 
with the estimated SFR, sSFR, and M*, suggesting that all 
star formation morphologies are present at all epochs irrespective 
of star formation rate and mass. 

• Simulated mergers among galaxies with similar masses and 
similar gas fractions typically predict similarly intense star forma- 
tion in both merger components, at odds with the observations: as 
opposed to 321^2^ of the observed mergers, as many as 59^6% 
of the simulated mergers show detectable star formation in both 
components. Together with the large fraction of early-stage, pre- 
coalescence mergers with star formation in just one component, 
this discrepancy supports the notion that z 1.5 mergers typically 
occur between galaxies with distinctly different properties such as 
gas fraction, mass and/or SFR. 
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APPENDIX A: PARAMETRISATION OF MERGER 
MORPHOLOGY (GM20CAS) 

To ease comparison between our visually selected merger candi- 
dates and the extensive literature selecting mergers based on empir- 
ically determined parametric classification schemes, we have esti- 
mated the Gini coefficient (G), the second-order moment of the 
brightest 20% of the fight (M20), the concentration (C), the asym- 
metry parameter (A), and the smoothness parameter (S) (e.g., Lotz 
et al. 2004, 2006, 2008b; Conselice 2003; Conselice et al. 2008, 
2009; Papovich et al 2005; Scarlata et al. 2007) for the objects pre- 
sented in this paper. We follow Lotz et al. (2008b) when estimating 
the individual parameters. 

In Figure Al the estimated GM20CAS parameters are shown 
for the 60 visually selected 3D-HST merger candidates (large black 
points) and the 292 visually discarded objects (small grey points). 
The histograms show the distributions of the individual parame- 
ters of the 60 merger candidates. The solid green and dashed red 
lines are from Lotz et al. (2008a,b) and Conselice et al. (2009), 



respectively. Parametrically selected mergers have been shown to 
preferentially lie above and to the right of these fines. We see that 
the majority of the visually selected 3D-HST merger candidates in- 
deed have A > 0.35 as expected. The G-M20 selection on the 
other hand seems to reject the majority of the merger candidates 
selected here. The selected merger candidates generally trace the 
parent population and does not distinct themselves clearly as would 
be expected for a 'clean' merger selection. This is most probably 
due to the complications that clumpy and irregular morphology of 
high- 2; objects and noisy images have on merger selections and in 
particular on 'blind' parametric classification schemes. 



APPENDIX B: THE 3D-HST EMISSION LINE (STAR 
FORMATION) MAPS 

In Figure B 1 we show the full sample of the 60 3D-HST emission 
line (star formation) maps. The objects have been sorted accord- 
ing to the estimated morphology of the star formation (SF type) 
described in Section 5. 
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